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Abstract 

We develop an athermal shear-transformation-zone (STZ) theory of plastic deformation in spa- 
tially inhomogeneous, amorphous solids. Our ultimate goal is to describe the dynamics of the 
boundaries of voids or cracks in such systems when they are subjected to remote, time-dependent 
tractions. The theory is illustrated here for the case of a circular hole in an infinite two-dimensional 
plate, a highly symmetric situation that allows us to solve much of the problem analytically. In 
spite of its special symmetry, this example contains many general features of systems in which 
stress is concentrated near free boundaries and deforms them irreversibly. We depart from conven- 
tional treatments of such problems in two ways. First, the STZ analysis allows us to keep track of 
spatially heterogeneous, internal state variables such as the effective disorder temperature, which 
determines plastic response to subsequent loading. Second, we subject the system to stress pulses 
of finite duration, and therefore are able to observe elasto-plastic response during both loading 
and unloading. We compute the final deformations and residual stresses produced by these stress 
pulses. Looking toward more general applications of these results, we examine the possibility of 
constructing a boundary-layer theory that might be useful in less symmetric situations. 
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I. INTRODUCTION 



The aim of this paper is to develop a theoretical description of the way in which the 
boundaries of voids in amorphous solids move under the influence of external loads. An 
obvious example is dynamic fracture. The evolution of cracks is traditionally described by 
linear elasticity. However, elasticity by itself does not provide a dynamical theory of crack 
motion but, instead, uses criteria such as energy balance to provide bounds on the velocities 
of straight cracks l|. This approach was challenged by one of us in More recently, 



S ta rtlllg with cloS e, y reiated idea, it „ shown inft that tahin g into accost the dynanhc 
degrees of freedom near a moving crack tip resolves various inconsistencies in traditional 
theories of linear elasticity fracture mechanics. Nevertheless, the theoretical approach of [2I 
is far from complete; it is still not an accurate and internally self-consistent description of 
free-boundary motion. 

We do not deal with the actual dynamics of cracks in this paper. Rather, we apply 
the shear-transformation-zone (STZ) theory developed in two recent papers 0, Q to the 
dynamics of a circular hole in an infinite medium. We have chosen this highly symmetric 
example in order to simplify the tensorial character of the theory, and to help in developing 
analytic methods that may remain applicable in less symmetric situations. As we shall see, 
the solutions of the circle problem provide a wealth of information about the ways in which 
stress concentrations near moving free boundaries induce plastic deformation. 

The STZ approach that we employ here is based on recent developments in the theory 
of dynamic elasto-plastic deformation of amorphous materials at low temperatures 4|. The 
theory was applied to a homogenous situation in |4| and the results were com pared to re- 
cent numerical simulations in la]. The STZ theory, originally proposed in 6|, |7|, deviates 
from conventional approaches [8( by focusing on the nature of the microscopic mechanisms 
for plastic deformation in amorphous media. These mechanisms are incorporated into the 
macroscopic description by internal state fields. An important observation is that the plastic 
strain itself cannot be one of these fields. For that to be true, the material would some- 
how have to encode information about its entire history of deformation starting from some 
reference state. This is not possible. More realistically, the memory of prior deformations, 
and the rate at which that memory is lost, must be encoded in the internal state variables 
and their equations of motion. The identification of the physically appropriate microscopic 
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variables is based here on the notion that STZ's are sparsely distributed, localized clusters 
of molecules that are especially susceptible to rearrangement in response to applied stresses. 
To take the STZ's and their role in the dynamics into account, one introduces a scalar field 
that describes the density of STZ's and a tensorial one that describes their orientation. 
Equations of motion for both the plastic strain rate and these internal state fields are ob- 
tained by assuming that STZ's change orientation (deform by finite amounts) at rates that 
depend on the stress. In addition, STZ's are created and annihilated at rates proportional 
to the rate at which the energy of plastic deformation is dissipated Q. 

An earlier STZ analysis of the circle problem Q| was based on a "quasilinear" version of 
the theory that used a conveniently truncated form of the transition rates 6|, |7| . This rela- 
tively simple theory captured linear visco-elasticity at small stresses, finite visco-plasticity at 
intermediate stresses, and unbounded plastic deformation at large stresses as a result of an 
exchange of dynamic stability between jammed (non-deforming) and unjammed (deforming) 
states of the system. Some memory effects were also described successfully. The quasilinear 
theory, however, had serious limitations; and therefore we base the present analysis on the 
recent version of STZ theory Q] that starts with more physically realistic assumptions. In 
4|, it was argued that the defining feature of a system at low temperatures is the constraint 
that, because thermal activation of STZ transitions is negligible or nonexistent, each molec- 
ular rearrangement occurs in response to an external driving force. No motion occurs in the 
absence of such a force, and no rearrangement moves in opposition to the direction of that 
force. In addition, this version of STZ theory includes the possibility that STZ's occur in 
many different types and with a range of different transition thresholds. 

We list here some important respects in which the present analysis differs from the pre- 



vious work on the circular hole problem in 9| and elsewhere 101 ]: 



Jl 

1. We employ here the fully nonlinear, athermal STZ theory described in [4j. 

2. The material time derivatives are expressed in a proper Eulerian formulation; and the 
theory is consistent with all relevant conservation laws and symmetries. 



3. The calculation described here takes into account the dynamics of t 



via an effective-temperature formulation, a feature that is shown in 
for self-consistency in these theories. 



lejiensity of STZ's 
to be essential 
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4. The focus here is on the transient growth regime that is relevant to dynamic situations 
such as fracture. Instead of considering only constant loads, we look at the elasto- 
plastic response to a stress pulse, and study the behavior during both loading and 
subsequent unloading. 

The organization of this paper is as follows: In Sec. [IT] we write the equations of motion 
and the boundary conditions for a circular hole in an infinite sheet of material subjected 
to a uniform radial stress at infinity. In Sec. IHH we review the elements of the STZ 
theory that are needed for this problem and specialize to cases of interest here. Sec. IIVI 
contains discussions of the specific situations to be considered in more detail, starting with 
the two different loading schemes to be used, and then the simplifying assumption of elastic 
incompressibility. Our numerical results for both constant applied stress and two different 
stress pulses are presented in Sec. [V] Sec. [Vljis devoted to the development of a boundary 
layer approximation and its comparison with the numerical results. Finally, in Sect. IVHt 
we offer a summary and further discussion. In Appendix [A], we show how to compute the 
threshold for unbounded growth of the circular hole under constant remote loading. 

II. EQUATIONS AND BOUNDARY CONDITIONS 
A. Strain, stress and rate of deformation 

We start by writing the full set of equations for a general two-dimensional elasto-plastic 
material. We define the total rate-of-deformation tensor 



where v(r,t) is the material velocity at the location r at time t. This type of Eulerian 
formulation has the advantage that it disposes of any reference state, allowing free discussion 
of small or large deformations. The price that we pay is that, in inhomogeneous situations, 
we need to employ the full material derivative for a tensor A, 




(2.1) 



VA 



— — h v ■ VA + A ■ uj — uj ■ A 

at 



(2.2) 



Vt 



where uj is the spin tensor 




(2.3) 
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For a scalar or vector quantity A the commutator with the spin tensor vanishes identically. 

Plasticity is introduced by assuming that the total rate-of-deformation tensor can be 
written as a sum of elastic and plastic contributions 

V<= el 

Dlf = _«+£$. (2.4) 
n X>t 13 y ' 

The linear elastic strain tensor efj is related to the stress tensor whose general form is 



CTij = -pSij + Sij , p = --Gkk • (2.5) 



The relation is 



where K and // are the two dimensional bulk and shear moduli respectively, is the 
deviatoric stress tensor, and p is the pressure. The equations of acceleration and continuity 
are 

Vv , . 

P Vt = V - (T = -Vp + V-s , (2.7) 

^ = - p V-v . (2.8) 

Consider now an infinite two-dimensional amorphous material with a circular hole cen- 
tered at the origin. The hole has radius R; and the system is loaded at infinity by a radial, 
uniform and possibly time-dependent stress a°°(t). For this configuration, the field of interest 
is the radial velocity v r (r, t), denoted here simply by v(r, t). The non- vanishing components 
of Dlf in Eq. (J2HJ are 

Dlf=^, = (2.9) 

or r 

Note that Dj° l satisfies the compatibility relation 

D rr = |: ( rD 0°o) ■ ( 2 - 10 ) 

The azimuthal velocity vanishes identically. The material time derivative in Eq. (12.21) be- 
comes 

We now turn to the equations of motion (12.71) . The symmetry of the problem implies 
that only the a rr and Oqq components of the stress tensor are non-vanishing. The tensor 
is diagonal, and thus we define the deviatoric stress s and the pressure p to be 

P = -7^{ a 8e + cr rr ) , s = sge = -s rr = -(agg - a rr ) . (2.12) 



Next we assume that the inertial term on the left-hand side of Eq. (12.71) is negligible and that 
density variations are also negligible; thus Eq. (12.81) is satisfied automatically, and Eqs. fl2.7p 



reduce to a single force balance equation of the form ll| 



^ + ^—^ = 0, (2.13) 
or r 

which, with Eqs.f l2.12]) . becomes 

I " 4I< A > . (2.14) 
The hole edge is a free boundary, i.e. it is traction-free; and therefore the boundary condi- 
tions are: 

a rr (R, t) = -p(R, t) - s(R, t) = 
p^-a°°(t), s^O as r-^oo, (2.15) 

where a°°(t) is assumed to vary slowly enough that the omission of the inertial terms in Eq. 
f l2~Tj) is still justified. 



B. Introducing plasticity 

Plasticity is introduced as in Eq. (l2.4p . According to Eq. (l2.6p . the components efj of the 
strain tensor satisfy 

,eZ P S 



2K 2fi 



ef g = -4- + — • (2.16) 
09 2K 2/i v ; 

The two-dimensional bulk modulus K is given by 

A- = + p , (2.17) 
1 - v* K ' 

where v* is the two-dimensional Poisson ratio. The plastic rate-of-deformation tensor Dfj is 
assumed to be traceless (incompressible plasticity), which in circular symmetry means 

Dg = -D£ = D pl . (2.18) 

The constitutive behavior is introduced by writing D pl as a function of the deviatoric stress 
s and the internal state fields. The material velocity at the edge of the hole v(R,t) is just 
the rate of change of the hole radius 

R(t)=v(R,t) (2.19) 



and therefore, using the second equation in (12. 9p . we see that the hole evolves according to 

(2.20) 



Once Dg°g(R, t) is specified, the last relation becomes an equation of motion for the boundary 
of the hole. 

C. The velocity field 



In order to explore the analytic structure of the problem, we use Eqs. (12. 4ft . (12.91) . (12.161) 
and (I2.18P to obtain 



v dv 
r dr 



1 

K 



dp 



P + v— 
dr 



s + v 



ds 
dr 



2D pl 



(2.21) 
(2.22) 



1 Vp 
~KVt 

v_dv = }_^s +2DP i = 1 
r dr fiVt it 

where dots denote partial time derivatives. It is useful next to eliminate s(r, t) and p(r,t) 

from these equations, and thus to obtain a differential equation for the velocity field v(r,t) 

that does not depend explicitly on those time derivatives. For that purpose, we operate 

with —(ii/Kr 2 ) (d/dr) [r 2 (-)] on Eq.(J2Z22D, differentiate Eq.flMU) with respect to r, add the 

results, and then use the partial time derivative of Eq. (12. 14H to obtain 

li \ ( d 2 v ldv v 



1 + 



K 



' d 2 v ldv 
dr 2 r dr 



1 

K 
// 



dv dp d 2 p Id 



r 2 v 



ds s 
dr 



AD pl dD pl \ 
K \ r dr I 



[2.23) 



This is a linear second order differential equation for v with coefficients that are functions of 
s(r,t) and p(r,t). This equation can be further simplified by defining v = v/r and operating 
with d r [v (•)] on Eq. (12. 141) . After simple manipulations we obtain 



K 



2s dv 
K dr 



K 



( ® 2y + %® v \ 
\ dr 2 dr J 

'AD pl dD pl \ 

r dr J 



[2.2A) 



This equation requires two boundary conditions. To obtain these, we note that the first 
boundary condition of Eq. (I2.15P implies 

V\p(R,t)+s(R,t)\ 



Vt 



. 



(2.25) 
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Then, using Eqs. (12.211) and (I2.22j) and simple manipulations, we obtain 

dv(R, t) 2Kv{R) - 2/2D pl (R, t) 



(2.26) 



dr (K + n)R 

To complete the derivation we need an expression for v(R) in terms of s and D pl . To that 
aim we integrate Eq. fl2.14l) from R(t) to oo and use the first boundary condition of Eq. (12.15j) 
to obtain the integral relation: 

<r°°(t) = 2 I" $ -^-dr . (2.27) 

jR(t) r 

Taking the time derivative of Eq. (12.271) . we obtain 

f°° Vs(r,t) dr = s(R,t)R a°°(fi 

jR{t) Vt r ' ' R 2 1 ' ' 

In order to use the last result we operate with 2 fii( 4 )(dr/r) on Eq. (12.22p to obtain 

1 r 00 Vs(r,t)dr r°° D pl 
fi JR(t) T>t r JR(t) r 

To evaluate v(oo), note that, for r 3> R, the stress field is purely elastic and therefore 

the solution of Eq. (12. 14ft is p = —a°° and s ~ r~ 2 . Then, the large-r behavior of Eqs. ( 12.211) - 

( T2T22|) predicts that 

v(oo) = —. (2.30) 
Substituting Eqs.(j2J28]) and f|273Q|) in Eq. fT2T29l we find 



1 VsirA) dr r°° D pl , . 



This equation plays a key role throughout the rest of this paper. Note that it already 
incorporates both the boundary conditions of Eq. (12.151) and that, together with Eq. (12.261) . 
it provides the boundary conditions that we need to solve Eq. (12.241) for the velocity field 
v(r, t). 



D. Putting it all together 

In order to put the equations in their final forms, up to specifying D pl , we define 

W(r,t) S *M. (2.32, 
8 



Then Eq. fl2.24p becomes 



1 + 



fi 



K 

2s{r,t) 



dr 



K 



a fAD pl (r,t) . ^dD pl (r,tY 
h s(r,t)- 



dr 



(2.33) 



which is to be solved using the boundary conditions given in Eqs. fl2.26p and (12.311) . The 
former becomes 



.R 



W{R,t) 



2K--2/iD pl (R,t) 



R 



(K + fi)R 

In terms of W(r,t), the equation of motion for the deviatoric stress s(r,t), Eq. fl2.22p . is 



2.34) 



- rW(r, t) 



s(r, t) + v(r, t) 



(9s (r, t) 
dr 



2D pl (r,t) . 



(2.35) 



The pressure p(r, t) can be calculated using Eq. ( 12. 14ft . According to Eq. ( 12. 32ft . the velocity 
field v (r, t) is given by 

Rr r 

v ( r , t) = — + r W(r', t)dr' . (2.36) 
R Jr 

The initial conditions are the linear elastic fields determined by the load applied at t = 



p(r,t = 0) = -<t°°(0) , 
8 (r,t = Q) = <t°°(0)^P- 



(2.37) 



III. PLASTIC DEFORMATION 



To complete this theory, we need to choose a specific plastic rate-of-deformation function 
D pl (s, ...) and equations of motion for the internal fields (denoted by the dots). As stated 
in the Introduction, we use the STZ theory described in reference [4|. Here, D pl (s, A, m) 
depends on the deviatoric stress s, on the normalized STZ density A, and on a field m that 
describes the local average STZ orientation. Although ordinary thermal fluctuations are 
absent at the low temperatures considered here, the concept of an effective disorder tem- 
perature remains essential 4j, [5]. The effective temperature T e ff characterizes the state of 
configurational disorder in the system. It equilibrates to the ambient temperature T at high 
T (relative to the glass temperature T g ), but may fall out of equilibrium at low T where 
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disorder is generated by the atomic-scale, configurational rearrangements that accompany 
mechanical deformations. We denote T e ff = (EsTz/ks) X-, where Estz is a characteris- 
tic STZ formation energy; then the STZ density is proportional to the Boltzmann factor 
exp (— Most importantly, the time variation of the STZ density is slaved to the dy- 
namics of x- 

The set of equations that describes plasticity in our theory {4] is: 

D pl {S,m,A) = -Aq(s,m) , (3.1) 

To 

pT = ^ g(a ' m) ( 1 -X e )' (3 ' 2) 

T) A 9 

_ s ~ g(s ~ m)(e -i/x_ A) ) (3.3) 

Dt T 

5t = ^-^ASg(s,m) (xoo - X) , (3-4) 
Pt c r 



where 



q(s,m)=C(s) f -— — m] , (3.5) 



and s = s/sj,, with s y being the dynamic yield stress [4j. Here Tq is an elementary time 
scale, e is a dimensionless constant of order unity, c is a specific heat expressed in units 
of k,B per atom and thus is of order unity, and Xoo is the asymptotic value approached by x 
during continuous deformation 4j. 

Note that many of the features of these equations are model independent. For example, 
Eqs. (13.11) and (13.51) identify m as a "back stress" that governs an exchange of dynamic 
stability at \s\ = 1 (i.e. when s = s y ). For \s\ < 1, the stable steady-state solutions of 
these equations occur at \m\ = 1 such that the system is jammed and D pl (s, m, A) = 0. For 
\s\ > 1, on the other hand, m = 1/s and there is a non- vanishing plastic flow (see [4j for 
more details). 

Essentially all of the specific material-dependent properties of this theory are contained 
in the function C(s). In the STZ model described in [4j, C(s) is related to the rate R(s) at 
which STZ's transform between their two orientations: 

s Mm + Btsi t m 



where R(s) is given by 



R(s) = - I (s-s a ) P(s a ;C) ds a , (3.7) 



Tq JO 
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Here, R(s) is an integral over a distribution P(s a ; () that reflects the fact that the STZ's 
can be of various types with different activation thresholds s a . The parameter ( that char- 
acterizes the distribution is the only material parameter here; it controls the width of the 
distribution. The mean value of the distribution is s y , the dynamic yield stress, that was 
shown in jj] to be the value of the deviatoric stress at which the system undergoes its dy- 
namic exchange of stability from jammed to flowing states. For finite values of ( there can 
be nonzero, sub- yield plastic deformation for \s\ < 1. This behavior is well documented in 
the literature 8J. We note that for s very small or very large, 

R(s) ~ S c+2 for 5-^0+ (3.9) 
R(s) -5-1 for 5 > 1 . (3.10) 

A basic assumption of the STZ theory is that the STZ's are sparsely distributed and 
only weakly interacting with each other. For this assumption to be valid, Xoo must small. 
Indeed, in independent applications of the theory to actual materials [yj] or simulations of 
such materials p, Q, Xoo was found to be of order 0.15 or less. For such values of Xoo, the 
density of STZ's, A = exp (—1/xoo), is of order 10 -3 , or very much smaller. We then notice 
that A appears as a rate-determining pre-factor on the right-hand sides of Eqs. (13. ip and 
(13. 4p . which govern the bulk system- wide variables D pl and \\ but A does not appear in a 
similar way in Eqs. (13. 2p or (13.31) . which pertain to the dynamics of individual STZ's. It 
follows that the plastic strain rate and the effective temperature respond much more slowly 
to changes in stress than do the internal fields m and A, and that the slow dynamics of 
the effective temperature controls the observable mechanical behavior of the system in most 
circumstances. We conclude that, as long as the characteristic time scale of the external 
loading is not significantly smaller than r exp (1/xoo), we can safely replace Eqs. ()3.2p and 
(13. 3p by their stationary solutions: 

(s/\s\ if|5|<l 
m = m (s) = < (3.11) 
1 1/s if \s\ > 1 

and 

A = e~ 1/x . (3.12) 
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The approximation in Eg. (13. lip tells us that \m\ = 1 for all \s\ < 1, i.e. the system 
"immediately" becomes jammed at small stresses. In fact, one can think of an experiment 
in which at time t = a virgin material with m = throughout is being used. Then 
in regions where s is small, m will remain very small for a long time, since C(s) is small. 
Notwithstanding, the precise value of m in regions of small s is not important for any of 
the results presented below. Therefore, we use Eqs. (13. lip and (13 . 12[) to reduce Eqs. (13. ip - 
O to 

D pl (~s } m,A) = ^e- 1 /*q (§) } q (s) = C(s) ( ± - m„(S) ) , (3.13) 

^ = —e- 1/ Xsq (s)( Xoo -x) ■ (3.14) 
Vt c r 

Once having made these simplifications, we also may go immediately to the limit £ — > oo 
and neglect any structure in the function C(s) for stresses less than or of order the yield 
stress. Without this approximation, our equations of motion would have predicted sub- 
yield deformations of magnitude proportional to C(s) exp(— 1/xoo), which we assume to be 
extremely small. Such deformations will not be negligible, of course, for soft and/or highly 
disordered systems, or for thermal systems at temperatures high enough to activate STZ 
transitions. But we do not need to consider such situations here, and therefore ( — > oo will 
suffice. In this limit, C(s) ~ 6>( | s | — 1) (\s\ — 1), where 9 is the Heaviside step function, and 

r , n/!1 ,x - l)(s sign(s) - 1) 
q (s) w Q(\s\ - 1) z (3.15) 

s 

The initial conditions in Eqs. (12.371) must be supplemented with an initial condition for 
the effective temperature: 

X (r,t = 0) = Xo, (3.16) 

where Xo describes a homogeneous state of disorder of the material before any load is applied. 
Because we are using Eq. (l3.1ip . there is no comparable initial condition for m. 

Eqs.flHl, UM-UM, QZM-UM, ASH-AM, and the initial conditions in 
Eqs. (12.371) and (13.161) define our problem. This system is specified by the following set of 
dimensionless parameters: eo, cq, Xoo, Xo, n/sy, K/s y and a given loading scheme 
a°°(t). The first three parameters characterize our plasticity theory; xo specifies the initial 
state of disorder of the material; and the last two parameters are the elastic moduli in units 
of the yield stress s y . Note that two parameters, r and s y , are scaled out as units of time 
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and stress respectively. The only length scale in the problem is the initial radius of the hole 
which, without loss of generality, we set to unity 



IV. LOADING SCHEMES AND THE INCOMPRESSIBLE LIMIT 
A. Loading schemes 

We now specify the loading scheme a°°(t). To make contact with earlier work we will 

look briefly in what follows at the case of increasing the load to a constant value, that is, 

f o"oi^— , for < t < 500r 
a QO (t) = \ 500r ° (4.1) 
I a , for t > 500r . 

Our primary interest here, as in 9j], will be to see how the region of plastic deformation 

forms, expands, and reaches a stable, stationary state for loading stresses <7o less than the 

threshold for unbounded growth of the hole (see Appendix lAl). 

We will pay greater attention to a class of situations in which the load takes the form 
of a stress pulse of finite duration, because we believe that such a pulse may be relevant to 
fracture dynamics. Plastic deformation of a material failing by crack propagation is localized 
(if it happens at all) near the fracture surfaces. At any given material point, deformation 
occurs only during a short time interval as the crack tip passes nearby. Specifically, a 
material element lying ahead of a crack tip experiences first an increasing stress as the crack 
tip opens, and then a decreasing stress as the load vanishes on the newly formed - and 
deformed - fracture surfaces. 

To simulate such a process in the circular geometry, we study a time dependent loading 
scheme in which the remote stress a°°(t) increases monotonically from zero to a peak stress 
a p > s y , and then returns monotonically to zero: 

f ia j^^ , 0<t<T, 
a°°(t) = { P T (4.2) 
[ 0, t > T. 

An example with T = 8000 r and o v = 2 s y is shown in Fig. [TJ 

Our choice of T in Eq. (14.21) is an order-of- magnitude estimate emerging from the analysis 
in Q]. For cracks whose speeds are governed by plastic dissipation, and when surface tension 
is negligible, that analysis implied that the crack tip blunts in such a way that the stress on 
its surface is always of order s y , and the tip radius is linearly proportional to the crack speed. 
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FIG. 1: A time dependent non-monotonic loading scheme corresponding to Eq. (|4.2p with T = 
8000to and a p = 2s y . 

Thus the time T that any material element spends in the fracture zone is roughly constant, 



of order the plastic relaxation time To exp ( 
estimate is consistent with the discussion in 



VXoo). If Xoo = 0.13, then T = 2200 r . This 
3j, where the plastic zone size was dynamically 



adjusted so that the time needed for the crack to pass this zone was of the order of the 
plastic relaxation time. Note that T is a very short time, of order picoseconds if To is an 
atomic time of order femtoseconds. That value of T is roughly an atomic spacing divided by 
a sound speed. These estimates lead us to believe that, although STZ plasticity is slow on 
atomic time scales, it is easily fast enough to be relevant to fracture. Clearly, this estimate 
of relevant time scales will have to be reexamined when the theory is applied to realistic 
models of fracture where surface tension is important. 



B. The incompressible limit 

It is convenient for both analytic and numerical simplicity to work in the limit of elastic 
(as well as plastic) incompressibility. Many materials are highly incompressible, as evidenced 
by large bulk moduli K or, equivalently, by the two-dimensional Poisson ratios v* being close 
to unity. Moreover, we have studied numerically the set of equations (12.141) . (I2.21l) - (l2.22p 
and (I2.32l) -( l3~8l) for various loading scenarios using 0.7 < v* < 1, and have found results 
that vary smoothly with v* without any noticeable singular behavior in the limit v* — > 1. 
Therefore, for analytic purposes, we can focus on the limit v* — >1 (K— >oo). 
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Note first that, in the incompressible limit, Eq. (I2.2ip reduces to 

v dv , 
- + 7T- = , 4.3 
r or 

which is in accord with the zero total compressibility trD tot = 0. Eq. (14.31) can immediately 
be integrated to yield 

, N R(t)R(t) 
v{r,t)= W W , (4.4) 
r 

where the boundary condition Eq. ( 12 . 1 9f) has been used. This explicit expression for the 
velocity field provides a major simplification because we no longer need to solve Eq. (12. 36ft 
numerically. 

Having the velocity field at hand, we can rewrite all the equations of the theory in a more 
explicit way. The sum of Eqs. (12.211) and d 2 . 2 2 [) becomes: 

RR^J_ds + J_RRd 1 + Dpl 
r 2 2/j dt 2ii r dr 

Eqs. (l3.13"]) - (l3.15p (with s = s/s y as before) are: 

=*.-«*»(.-), qo( g) =e(|i|-l) (4 . 6) 

T S 

dx RRdx 2e _ x/ 

^7 + -z- = e /x sq {s) (Xoo ~ X) ■ O 4 - 7 

at r or cq r 

The derived relation, Eq. (l2.3ip . becomes: 



roo dp 1 , &°°(t) 
2 / dr + 



R _ jr r 2/i 
R~ 1 s(R,t) 



(4.8) 



Note that the pressure p(r, t) does not appear in the final equations for the incompressible 
limit, but is computable as before from the deviatoric stress field s(r,t) using Eq. (12.141) . 



V. RESULTS 

We are now ready to study numerically Eqs. (|4.5j) - (l4.8j) in detail. We consider both the 
constant load and the time dependent loading scenarios, Eqs. (14. II) and (14.21) respectively. 
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FIG. 2: The r dependence of the dynamical variables at different times in the constant load case. 
The radius of the circle is shown in the upper left panel without the initial elastic expansion. 
At longer times the radius saturates. Note the relative small increase in radius. The stress is 
shown in the upper right panel, and the fields D pl (r,t) and A = e~ l ^ x in the two lower panels 
respectively. The parameters chosen are eo = 1, Co = 1, Xoo = 0.13, xo = 0.1, [i = 50s y 
and a 00 = 2s y . The solid line corresponds to t = 2000to, the dashed line to t = 4000to, the dotted 
line to t = 6000ro and the dotted-dashed line to t = 8000to. 

A. Increase toward a constant load 

The results of numerical simulations for the case of an increase toward a constant load, 
&o/s y = 2, are shown in Fig. [21 In this loading scheme, we know that there is a maximum 
value of <7°°, say a th , above which the hole expands indefinitely. (See Q| and earlier references 
cited there.) For completeness, in Appendix El we show how to compute a th in the current 
version of STZ theory. The result, for the parameters used here, is a th « 5^, which is 
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substantially larger than the value of <To chosen for this illustration. 

Note that the radius of the hole, R(t), first increases elastically to Rq = exp(<7o/2yu) = 
1.02 (not shown in Fig. [2]), then grows (by less than one percent in this example), and tends 
toward a constant value at large times. The stress s is proportional to the elastic solution, 
1/r 2 , for r outside the plastically deformed region; but, in accord with conventional plasticity 
theories, s decreases in time toward the yield stress s y inside that region. The 
eventually extends out to a radius R\. The usual quasistatic estimate of R\ [8|, 
by assuming s{r) = s y for R < r < R x , and s(r) = s y (Ri/r) 2 for r > R x . Inserting this 
s(r) into the integral on the right-hand side of Eq. (12.271) and setting a°° = <jq = 2 s y on the 
left-hand side, we find R x /R = exp[(l/2) (cr /s y — 1)] = 1.65. Thus, even in cases where 
the plastic deformation is small in comparison to the elastic displacement, the plastically 
deformed region may be quite extensive. Finally, note that the STZ density A = e~ 1//x 
becomes substantially larger than its initial value throughout this plastic region. 



astic region 
3 is m ade 



B. Time dependent non-monotonic loading 

In Figs. [3] and H] we present numerical results for the case of time-dependent non- 
monotonic loading as shown in Eq. (14.21) and Fig. [TJ Here we have chosen T = 8000 r 
and two different values of the peak stress, <J p /s y = 2.0 and 4.0. The first of these is small 
enough that the radius of the hole changes at its maximum by only about two percent, 
consistent with the effect of constant loading at the same stress shown in Fig. [2j Only a 
small part of that transient change in the radius is irreversible plastic deformation; most of 
what we are seeing is elastic expansion and contraction of the system as a whole. However, 
this relatively moderate scenario may be the more realistic of the two sets of time depen- 
dent simulations for describing failure in strong structural materials where, in keeping with 
conventional analyses, we do not expect stresses to be much larger than the yield stress. 

The bigger peak stress causes substantially more plastic deformation, as seen in the 
comparatively large, irreversible change in the hole radius that remains after the system 
has been fully unloaded. This peak stress may be unrealistically large; but showing it here 
illustrates some features of the plastic response more clearly than the small-stress case. Note 
that the large applied stress induces a complex unloading sequence in which the stress near 
the hole becomes so large and negative that it drives strain recovery via reverse plastic 
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FIG. 3: The time dependent non-monotonic load case. All the material parameters are the same 
as those for Fig. [21 while the load is given by Eq. (|4.2j) with T = 8000to and a p = 2s y . The solid 
line corresponds to t = 2000to, the dashed line to t = 4000ro the dotted line to t = 6000to and the 
dotted-dashed line to t = 8000ro. 

deformation. 

The dynamic response of this system, in both examples, is characterized by the emergence 
of a residual stress in the plastic zone near the hole edge. Note that s exhibits the sharp 
transition between plastic and elastic regions that we expect as a result of the special form 
of the constitutive law chosen in Eq. fl3.15p . and that the STZ density A becomes large in the 
plastic region. The qualitative picture that we take from these results is that of a plastically 
deformed region near the hole edge (or near any defect boundary in the general case), outside 
of which s retains its elastic solution proportional to 1/r 2 , but within which s decreases and 
may change sign, consistent with an irreversible outward displacement of the material near 
the hole. 
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FIG. 4: Same as preceding figure but with a p = 4:S y . 
VI. BOUNDARY-LAYER APPROXIMATION 

A principal objective of this investigation has been to study the dynamics of time- 
dependent plasticity in the neighborhood of a symmetrically loaded circular hole, and, from 
this study, to learn how to construct realistic approximations for less symmetric situations. 
We have seen in the preceding Section that, under certain circumstances, plastic deforma- 
tion remains localized in a relatively narrow zone near the hole. This observation leads us to 
believe that some sort of boundary- layer approximation might have at least limited validity 
in interesting applications such as asymmetric void growth or fracture. In what follows, we 
explore one class of such approximations. 

The numerical examples illustrated in Figs. [3] and H] have several common features that 
need to be captured in any useful approximation. Most importantly, the dynamic behavior 
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is qualitatively different in what we call the "active plastic phase" - roughly speaking, the 
time interval in which the hole is growing by plastic deformation - than it is in the "elastic 
unloading phase" where plastic deformation ceases. During the active plastic phase, as soon 
as cr°° exceeds s y , plastic deformation occurs in an active plastic zone, R(t) < r < R\(t), 
where Ri(t) is the radius at which 

s(i2i(*),*) = s„. (6.1) 

During this phase, no irreversible deformation has occurred for r > Ri(t), and thus the 
stress field in this outer region is simply 

s(r,*) = V^T^ P(r,t) = -o°°(t), for r > R x . (6.2) 

This behavior changes as elastic unloading begins. As seen in Figs. [3] and HI R(t) reaches 
its peak considerably later than when the load cr°°(t) peaks at t p = 4000 r . Between t p 
and the time at which R(t) reaches its maximum, and for some time after that, the hole 
continues to grow plastically while the system as a whole contracts elastically. During that 
interval, the stress at the boundary of the hole, s(R, t), drops rapidly toward the yield stress 
s y ; while s(Ri,t), by the definition of Ri, remains equal to s y . The transition between the 
active plastic and the elastic unloading phases that is important for our purposes occurs 
when s(R,t) passes downward through s y at a time that we call t\. 

To see what happens near ti, we show in Fig. [51 for the case a p = 2s y , a sequence of 
graphs of s(r,t) and D pl (r,t) as functions of r/R(t) for four equally spaced times starting 
at t p = 4000 To, and ending at t — 5500 r (just later than t\ = 5200 r ). Note that at 
t = 5500 ro D pl has already vanished. The crucial observation is that D pl (r,t) extends 
across the whole active plastic zone, R(t) < r < Ri(t), until t = t\, at which time, at least 
to a rough first approximation, it vanishes almost uniformly everywhere. 

The elastic unloading phase, starting at about time t\, evolves as an entirely elastic 
response to the decreasing driving force, conditioned by the material displacements that 
occurred during the plastic loading phase. During this part of the unloading process, there is 
no longer an active plastic zone, but there remains a plastically deformed region - a "process 
zone" - within which irreversible material displacements give rise to residual stresses. The 
process zone extends from R(t) out to the most distant material point at which plastic 
deformation occurred, that is, out to advected back toward R(t) by elastic relaxation. 
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FIG. 5: The profiles of s(r,t) and D pl (r,t) after the external stressthe same as those for Fig. [3l 
The solid line corresponds to t = 4000to, the dashed line to t = 4500to the dotted line to t = 5000ro 
and the dotted-dashed line to t = 5500to. 

Elastic unloading persists at least until t = T when the remote load cr°°(t) vanishes. As 
seen in Fig. HJ the stress near the hole may become so large and negative that it drives 
plastic strain recovery. This behavior occurs only for large peak stresses, and does not occur 
when the remote stresses are smaller than about 2.5 s y in the present model. We will not 
include plastic strain recovery in our boundary- layer analysis, and simply will point out 
where it is missing. 

Our strategy for developing a boundary-layer approximation is to write equations of 
motion for quantities defined only near the boundary, i.e. within the process zone, and 
then to deduce the boundary motion from these local quantities instead of solving for fields 
everywhere in the system. In our case, the relevant boundary is just the hole radius R(t). 
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The rate of change of any quantity A(R, t), defined on r = R(t), is 

dA(R,t) dA(R,t) b dA(R,t) 

~^ = ^^ + R ^—- (6 ' 3) 

Since R = v(R,t), the time derivative in Eq. (16.31) is just the material time derivative 
defined in Eq. (12.111) . The equations for s(R,t) and A(R,t), i.e. Eqs. ( 14.51 - H~Tj) evaluated 
at r = R(t), are: 

D P\R,t) = ^e- 1 M R ^q (R,t); q (R,t) = q (s(RM (6-5) 

TO 

= ^ e"^^) s(R, t) q (R, t) [ Xoo - X (R, t)} . (6.6) 
dt c r 

We need one more relation to close this set of equations. The obvious candidate is 
Eq. (14.81) . which already involves only quantities defined in the process zone. This equation 
demands special attention, however, because it has been derived using exact mathematical 
relationships that are unique to circular symmetry, and because it explicitly involves the 
remote driving force a°°(t) which entered the analysis via those exact relationships. We 
need some such relation to determine the coupling between the remote driving force and 
the boundary. If the boundary- layer strategy is to be successful, we ultimately will have to 
interpret a°°(t) as a quantity emerging from solutions of a more general elasticity problem. 
For present purposes, however, we accept Eq.1 14.8p as written, except that we simplify it by 
neglecting the term s(R,t)/fi ~ s y /fi <C 1 in the denominator on the right-hand side. Thus 
we write: 

R rRi D pl (r,t) , &°°(t) , , 

^ ° 1 ' dr + ^ . (6.7) 



R Jr r 2/i 
Note that the original integral in Eq. (14. 8p involves the plastic rate of deformation D pl (r,t) 
which vanishes for r > R\. Therefore, in Eq. (16.71) . we have explicitly inserted Ri as the 
upper limit of integration. 



Next we use Eq. (16.71) to obtain an approximate expression for R/R, and then use that 
relation to eliminate R/R on the left-hand side of Eq. (16.41) in favor of s(R, t). To do this, we 
need an approximation for the integral over D pl (r, t) that appears in Eq. (l6.7p . The numerical 
results shown in Fig. [5] suggest that, throughout the active plastic phase, D pl (r,t) varies 
almost linearly, from its value D pl (R,t) at r = R(t) to zero (by definition) at r = Ri(t). 
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Thus we write 

To estimate R±(t) we invoke Eg. (12.271) . which is another exact relationship between the 
circularly symmetric stress field s(r, t) and the remote driving force a°°(t), and is subject 
to the same concerns that we expressed above about our use of Eq. (I4.8j) . So long as we 
remain in the active plastic phase, and can assume that Ri(t) remains the outer boundary 
of the material region in which D pl (r, t) is nonzero, we can use Eq. (16.21) for the stress outside 
r = Ri(t), and write Eq. (12.271) in the form 

= 2 J m '^f±dr + Sy . (6.9) 

Then, in the spirit of our approximation for D pl (r,t) in the active plastic region, we make 
a linear approximation for s(r,t), using the fact that this stress is equal to s y at Ri, and 
allowing s(R,t) = s(R(t),t) to be an as yet undetermined function of time. Specifically, for 
R(t) < r < Ri(t) and for times t <t\ such that s(R,t) > s y , we have: 

R, (t) — r 

s(r,t)^s y + R £ J _ (s(R,t)-s y ) , s(R,t)>s y . (6.10) 
Inserting Eq. (16.101) into Eq. (16.91) . and linearizing in (i?i — R)/R, we find, again for t <t\, 

(If- 1 ) p-") 

Note that when cr°°(t) is increasing and is smaller than s y , there is no active plastic zone at 
all, i.e. Riit) = R{t). Using Eq.flES), Eq.(EZD becomes 

We emphasize that, in this approximation, the end of the active plastic phase occurs at 
the time i x when s(R,t 1 ) = s y . By definition, D pl (R 1 ti) = 0, and therefore D pl (r,ti) = 
throughout R(ti) < r < Ri{t\). Apart from the possibility of reverse plasticity, D pl (r,t) 
vanishes at all later times, t > t\. It follows that Eq. (16.12j) . with a properly interpreted 
D pl (r,t), is a valid approximation at all times, during loading and unloading, because - 
apart from the possibility of reverse plasticity - D pl (R,t) vanishes whenever s(R,t) < s y . 

Our boundary- layer calculation is straightforward from here on. We insert Eq. (16.121) on 
the left-hand side of Eq. (16.41) . thus obtaining a nonlinear, first-order differential equation 
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FIG. 6: A comparison between the exact solution (solid line) and the prediction of the boundary 
layer theory (open circles) for a p = 2s y . 

for s(R,t), which can be solved using Eqs. fl6.5H - fl6.6j) . Our results are shown in Figs. [MEl 
where we compare the predictions of the boundary layer theory for R(t), s(R,t), D pl (R,t) 
and A(R,t) with the exact solutions for cr p /s y = 2, 3 and 4. The comparison is excellent 
for Up/sy = 2, even though Ri(t)/R(t) — 1 is about 0.5, which is not that small. The results 
are acceptable also for cr p /s y = 3, but deviations grow rapidly for higher values of u p /s y , as 
expected. Note that in Figs. [7] and [8] the effect of plastic strain recovery is noticeable for the 
exact solution but has not been included in the boundary- layer approximation. Nevertheless, 
the values s(R(t)) < —s y and D pl (R(t)) < are well approximated in this period of plastic 
strain recovery. 

We can push the comparison between the boundary- layer approximation and the exact 
results further by computing the residual stresses in the process zone. We do this as follows. 
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FIG. 7: A comparison between the exact solution (solid line) and the prediction of the boundary 
layer theory (open circles) for a p = 3s y . 

We note first that Eg. (I6.4j) can be used for any radius R'(t), not just the radius of the hole 
R(t). At times later than t±, the quantity D pl (R', t) appearing here vanishes for any R' . We 
also can compute the quantity Rl / R' on the left-hand side of Eq. (16. 4p using our knowledge 
of R(t)/R(t) and the fact that the material in the region between R and R 1 is incompressible 
(R' R' = RR). Finally, we know that s(R',ti) = s y for all R' between R and R\\ so we can 
use Eq. (16.41) to compute the advected value of s(R'(t),tj at any later time t, and then plot 
this value as a function of the advected position r = R'(t). The functions s(r, t) computed 
in this way, for <7 p /s y = 2, are shown in Fig. [91 Outside the process zone, we have used the 
elastic solution with s oc 1/r 2 . The agreement between the full solution and the boundary 
layer predictions lends further support to the philosophy presented in this section. 
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FIG. 8: A comparison between the exact solution (solid line) and the prediction of the boundary 
layer theory (open circles) for a p = As y . 

VII. DISCUSSION AND SUMMARY 

We have presented a derivation of the equations of motion for a moving free boundary, 
in this case a circle, using linear elasticity and the athermal STZ theory of amorphous 
plasticity. These equations were solved numerically and yielded some interesting results. We 
looked especially at loading scenarios in which the material was subjected, at large distances 
from the hole, to stress pulses whose durations were comparable to the plastic relaxation 
time for the STZ mechanism. For the range of material parameters chosen here, which 
we believe to be characteristic of realistic amorphous solids, the irreversible displacement 
of the boundary of the hole was small. However, the width of the region in which plastic 
deformation occurred, as evidenced by increased internal disorder and residual stresses, 
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FIG. 9: Comparison between the full solution (solid line) for s(r,t) and the prediction of the 
boundary layer approximation (dashed line) for a v = 2 s y and for various values of t. 

became comparable to the radius of the circle for stress pulses whose peak strengths a p were 
twice the yield stress s y . We also found that strong enough stress pulses, o v > 2.5 s y induce 
reverse, plastic strain recovery near the edge of the hole. 

The theory as written is rather cumbersome even for the circular symmetry considered 
here. It is therefore with some relief that we have found that a boundary layer approxima- 
tion, using only quantities defined on the circumference of the circle and their coupling to 
the remote stress, succeeds in capturing the plastic effects quite well. Specifically, despite 
several apparently serious oversimplifications, this approximation accurately reproduces the 
irreversible displacement of the radius of the circle and the residual stress found in its neigh- 
borhood. We propose that such boundary layer theories might find useful applications in 
less symmetric situations, not the least interesting being fracture dynamics. 
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A next issue on our agenda is the dynamic stability of our circular solutions. It will be 
interesting to learn whether our equations of motion predict that the growing circular hole 
becomes unstable against symmetry-breaking perturbations, for example, whether it forms 
fingers that might evolve into cracks, or whether a perturbed circular solution can survive. 
The transition between these two possibilities, if found, might shed light on the brittle-to- 
ductile transition in visco-plastic materials. Along the same lines, it will be important to 
learn how strongly the circular behavior itself, and the stability of circular symmetry, depend 
on specific details of our model such as elastic compressibility, initial states of disorder, or 
the absence of thermal effects. 
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APPENDIX A: UNBOUNDED GROWTH 

One characteristic stress that emerges in the circular-hole problem is the threshold a th 
above which the hole grows without bound under conditions of constant loading. In this 
section we explain how this threshold is estimated in the athermal STZ theory. 

Because the only length scale in this problem is the radius of the hole, R(t), we expect 
that all spatially varying (r-dependent) quantities occurring in a uniformly (exponentially) 
expanding solution of our equations of motion will be self-similar functions only of the ratio 
r/R(t). For convenience, we choose the variable £ = R(t)/r. These self-similar solutions 
are characterized by R/R — u, where u > depends on the loading and vanishes at the 
threshold o th . 

In fact, these solutions invalidate the assumptions of our theory for large times since 
the exponential growth of the hole implied by R/R = uj will inevitably lead to very large 
velocities at large r, in contradiction with our omission of the inertial term in Eq. (12.71) . 
When this happens in the full, inertial theory, the self-similar solutions break down with 
the appearance of perturbations propagating at a characteristic wave velocity. Nevertheless, 
as we are interested in obtaining an estimate for a th , the self-similar approximate solutions 
for finite times are still useful so long as we understand that we are actually dealing with a 
system in which the remote tractions are being applied at a large but finite distance from 
the hole. 

For present purposes, we solve our equations of motion for a self-similar deviatoric stress 
s(£;u;) in the limit t lj — ► 0. Once s(£;a>) is known, we can use Eq. (I2.27P to obtain the 
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following expression for a 



th 



^c2f 1S ^dt for (Al) 

Sy JOt, 



Here we assume that the exponential growth at a small u allows the system to approach the 
self-similar solution u) in times when the velocities are still smaller than the typical wave 
velocity s(£;a;) can be computed using Eqs. ( 14. 51) - (14. 71) . First we note that the material 
time derivative of Eq. (12.111) translates to 



2. = u £h-?)9. (A2) 



Then Eqs. (gSD-gZD read 



d s 

t u£(1-£ 2 )— = 2/ir u;£ 2 - 2/2e e' 1/x q(s } m;C) 
din 

r ^£(l-£ 2 )— = 2g(S,m;C) (1-sm) 

r,uji{l-e)% = — e- 1/x 5g(5,m;C)(Xoo-x) • 
Ok, c r 



(A3) 



Also recall that R/R = to. The initial conditions are 

5(0) = 0, m(0) = 0, x(0) = Xo- (A4) 

We have integrated these equations numerically in the limit t q uj — > with e = 1, c = 
1) Xoo = 0.13, Xo — 0.1, /V s ?/ = 50 and for various values of £ < 10. Then we have used 
Eq. (]Aip to obtain cr' h « 5 s^, and have checked that this result is only weakly dependent 
on C- 



30 



